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■ Abstract 

(N 

Using the formalism of extended = 4 supersymmetric quantum mechanics we 
consider the procedure of the construction of multi-well potentials. We demostrate 
the form-invariance of Hamiltonians entering the supermultiplet, using the pre- 
sented relation for integrals, which contain fundamental solutions. The possibility 
(— I ' of partial = 4 supersymmetry breaking is determined. We also obtain exact forms 

of multi-well potentials, both symmetric and asymmetric, using the Hamiltonian 
I of harmonic oscillator as initial. The modification of the shape of potentials due 

J> ' to variation of parameters is also discussed, as well as application of the obtained 

^ ■ results to the study of tunneling processes. 

\o 1 Introduction 
O 

^ ' Hamiltonians with multi-well potentials are of interest for classical, as well as quantum 
systems. Among the most interesting applications of multi-well potentials in classical 
dynamics is the problem of transitions of particles from one local minimum to another 
under the influence of different types of noise. This problem is common for many physical 
5^ ! systems from Universe to microcosm. History of calculations of transition rates from local 
minimum to another take rise in the famous paper of Kramers [I] and after 70 years is still 
far from completion. Especially interesting is observed effect of amplification of transitions 
over the barrier under the action of weak time dependent periodic signal, which is called 
stochastic resonance [21 13] . 

Dynamics in multi-well potentials in quantum mechanical case is not less interesting. 
That's enough to mention the tunneling effect, which is especially interesting in the con- 
nection with wide-ranging research of trapping of atoms of alkali metals in superfiuid state. 
Moreover, when more than one barrier exists, the resonant amplification of tunneling rate 
is possible, which is known as resonant tunneling [3]. Experimental observation of this 
phenomenon in superconductive heterostructures |5] serves as a basis for construction of 
resonant tunneling diode. It is important to have mechanism to change the parameters of 
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potential (locations of minima, heights of a barriers), to determine the conditions of this 
phenomenon. 

Above mentioned systems, classical as well as quantum, are joined by methods of 
theoretical analysis. Fokker-Planck equation (FP) is one of the most important kinetic 
equations, that describe the dynamics of stochastic systems [6l |7]. One of the useful 
methods for solving FP equation is the eigenfunction expansion method. Due to formal 
similarity of FP and Shrodinger equations, an eigenfunction expansion method, analo- 
gous to a bound state expansion treatment of the Schrodinger equation, is very useful for 
the Fokker-Planck equation. In particular, in the case of bistable stochastic system, the 
knowledge of full set of wave functions and eigenvalues of corresponding quantum Hamil- 
tonian completely determines time evolution of the solutions of FP equation. Moreover, 
this allows to make a conclusions about the dynamics of corresponding metastable sys- 
tem. Theoretical analysis of the processes in multi-well potentials is complicated due to 
the fact that existing models operate usually with piecewise potentials (e.g constructed 
from rectangular or parabolic wells and barriers), which are far from real potentials. Aside 
from this, the wave functions and spectrum are unknown in such potentials, which implies 
the numerical analysis of their properties. 

Exact solutions of FP equation in this approach are closely connected with existence of 
exactly solvable quantum-mechanic problems, which number increases significantly during 
last years [HI [9] . It is important to mention that the first meaningful example of construc- 
tion of the significantly nonlinear models of diffusion in bistable system is presented in 
papers pU| [TT] . The formalism of Darboux transformation, used in this papers, have be- 
came the one of the most widely used methods of construction of isospectral Hamiltonians 
in supersymmetric quantum mechanics(SUSY QM). Existence of exactly solvable [12] and 
partially solvable [13J quantum mechanical models with multi-well potentials must fill the 
gap in theoretical analysis of above mentioned processes and serve as a more realistic 
approximation in their research. 

The aim of the present paper, first of all, is to determine the relations between dif- 
ferent types of potentials using the framework of the extended supersymmetric quantum 
mechanics (A^ = 4 SUSY QM) [TH [15j. Also, we research the possibility of changing 
parameters of potentials in wide limits and obtaining the exact expressions for spectrum 
and wave functions. This is especially important when potential has more than one local 
minimum and, thus, "resonant" tunneling is possible. Obtained expressions would be 
used for derivation of new exactly solvable stochastic models, which is natural extension 
of previous research in single- well potentials [16] . 

In the first chapter we give short explanation of A^ = 4 SUSY QM structure and 
procedure of construction of isospectral Hamiltonians with additional states above the 
ground state of the initial Hamiltonian. We especially emphasis the emergence of multi- 
well potentials, based on the general properties of the solutions of auxiliary equation of 
Shrodinger type. In particular, we discuss the possibility of partial supersymmetry break- 
ing in A^ = 4 SUSY QM. In the next chapter we give the expressions, that allows to 
analyze the obtained potentials and wave functions, regardless of concrete form of initial 
Hamiltonian. In particular, this allows to generalize the concept of form-invariant poten- 
tials [I?] and calculate in general normalization constants of zero-modes wave functions. 
Using the harmonic oscillator Hamiltonian as initial, we obtain the expressions of the po- 
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tentials of isospectral Hamiltonians and wave functions for the case of exact and partially 
broken supersymmetry. The analysis of conditions for emergence of double and triple- 
well potentials and ways for varying of their form in wide range is provided in the third 
chapter. In conclusion we summarize the main results and their possible applications. 

2 Isospectral Hamiltonians of = 4 SU SY QM with 
additional states 



Extended = 4 SUSY QM [HI |T5] is equivalent of second-order polynomial SUSY QM 
(reducible case) [TSl dS] and assumes the existence of complex operators of sypersymme- 
tries Qi{Qi)aRd Q2{Q2), through which the Hamiltonians H^^ could be expressed. 
Hamiltonian of = 4 SUSY QM has a form(^ = m = 1): 



CTl 



a 



(1) 



Vi{x)) 



Vi{x) 



W\x) + -af 



l^ip' + V^x] 
1 l^^W"{x) 



W'ix) 



where W{x) is a superpotential and a^^^ - matrices, which commute with each other and 
have eigenvalues ±1, = as ® 1, uf^ = 1 ^ as. 

Supercharges Qi of extended supersymmetric quantum mechanics form an algebra: 



{Qi, Qk} = 26ikH, {Qi, Qk} = {Qi, Qk} = 0, t,k = l,2 
Qi = a^\p + ^V,+^{x)), Qi = a? {p - iVi+,{x)) 



(2) 



where Vs{x) = Vi{x), a^^'^ = a± 



1, -i^^ 



1 ® a±. 



Hamiltonian and supercharges act on four-dimensional internal space and Hamiltonian 
is diagonal on vectors ip'^^{x,E), where ai, (T2 - eigenvalues of (T3 , Supercharges 
QiiQi) act as lowering (raising) operators for indexes ai, a^- It is convenient to represent 
the structure of Hamiltonian and connection between wave functions in diagram form: 



Qi 




Obviously, due to commutativity of operators Qi and Qi with Hamiltonian, all ip'^^ {x, E) 
are eigenfunctions of Hamiltonian with the same eigenvalue E. The only exception is the 
case, when wave functions turn to under the action of generators of supersymmetry. 
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Construction of isospectral Hamiltonians in the framework of = 4 SUSY QM is 
based on the fact that four Hamiltonians are combined into supermultiplet H^^. Never- 
theless, it should be noted that due to symmetry o"i ^ (J2 — H~^^ = HZ%^ only three 
of them are nontrivial. Let's note, that such relation in the case of higher- derivative 
SUSY [19] establishes the correspondence between quasi-Hamiltonian and operators of 
Shrodinger type and is identical for any superpotentials. The procedure of construction 
of isospectral Hamiltonians, when ground state is removed from the initial Hamiltonian, 
is treated in [20] in detail. We will consider the construction of isospectral Hamiltonians 
by addition of states above the ground state of initial Hamiltonian. Similar procedure 
was already performed (e.g. in [20]), but the distinctive feature of present research is the 
obtaining of general results without specification of concrete form of initial Hamiltonian. 

Let's consider the auxiliary equation: 

Hip{x) = eip{x) (3) 
As initial let's take one of the Hamiltonians 

^"f = lip - ^(riV2ix)){p + t(T^V2{x)) +e = 
= lip - ia2Vi{x)){p + ia2Vi{x)) + e 

where e - the so-called factorization energy. Hereinafter the energy is measured from 
e. Strictly speaking, the supersymmetry relations with supercharges of the form ([2]) and 
expressed through superpotential W{x) are satisfied not by Hamiltonian H itself, but 
by a shifted by an amount of e Hamiltonian H — e. However, due to commutativity of 
supercharges with constant e, the relations between wave functions of H and H — e are the 
same. When fixing operator H^^, the form of W{x) depends on the choice of factorization 
energy e, thereby the Hamiltonians H'^-^^^, H^^'^, HZ^f also have nontrivial dependence 
on e. 

When e < Eq (where Eq is the ground state energy of initial Hamiltonian), the aux- 
iliary equation ([3]) has two linear independent solutions (pi{x,e), i = 1,2, which are 
nonnegative and have the following asymptotic behaviour [2T|: when x — > —00 ^pi{x) — )■ 
+00 {ip2{x) —J- 0), and when x — t- +00 v'i(x) — > {ip2{x) — j- +00), i.e. with appropriately 
chosen constants the general solution has the form ip{x, e, c) = N{{pi{x, e) + ap2{x, e)) {N 

- normalization constant) and has no zeros on all axis. Thus, the function ip{x,e,c) = 
AT-i 

— is finite and could be normalized at every concrete choice of e and c. Let's note 

ip{x,e,c) 

that with concrete values of e and c ip{x,e,c) could have local extrema. In this case the 
natural choice of initial Hamiltonian is or (which are identical due to symmetry 
of H^^ under (Xi ^ 0-2) • Superpotential in this case has the form: 

1 ^ 

W{x, e, A) = -- ln(l + \J dt^'^{t, e, c)) (5) 

where A, Xj - two new arbitrary parameters, but one of them is inessential, because it 
gives an additional contribution to W{x). All Hamiltonians that form a supermultiplet 
have nontrivial dependence on these parameters. 
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To consider the connection between Hamiltonians from supermultiplet, let's take 
as initial. Denoting as ip^{x,E) the solution of equation 



H-i,-{x,E) = Ei,-{x,E) 



(6) 



and using the first representation of Hamiltonian we obtain the following relation 
between, ipz{x^E) and initial expressions: 



H 



ip_{x,Ei 



1 



dx'^ 



\n(p{x,e,c), 
W [ip-{x,Ei),(p{x,e,c)] 



2{E, 



ipZix,E = 0) 



N 



-1 



^{x,e,c) 



(p{x,e,c) 



(7) 



The new state with E = (with energies measured from e) appears in Hamiltonian, which 
by definition has normalized wave function. Normalization of wave functions of excited 
states is preserved (for discrete spectrum). 

Using the second representation H^^ = -{p — i<72Vai{x)){p+ia2Vo-j^{x))+e and identity 

= the relation between H^, ip^lx, E) and initial Hamiltonian could be obtained: 



Hl = H 



dx"^ 

%l)X{x,E = 0) 



In 



(p{x,e,c) 



\ 



ipl{x,Ei 



2{Ei - e) 



— + — In 

doc ditjc 



l + XJ dt <^2(t,e,c) 
{l + XJdtip^it,e,c)) 

Xi 

d , ip{x,e,c) 



\ 



\jdt^\t,e,c)) 



tp^{x,Ei) 



It is important to mention that normalization of wave function iIj'X{x,E = 0), as in the 
case of one-well potentials, could always be performed at any (p{x,e,c) and A using the 
following expression in normalization condition: 



N^'ip\x,e,c) 



(l + AAr-2 / dt^^{t,e,c)y 



dx 



(l + AAr-2 J dt^^{t,e,c)) 



which easily gives the relation between normalization constants: A^^^^ = (1 + A)A^~^. 
Normalization of ip1^{x,Ei) is the same as for tlj^{x,Ei) at any choice of ip{x,e,c). The 
usage of superpotential with ip{x,e,c) corresponds to exact supersymmetry, which 
leads to existence of zero-modes in Hz and H^. Existence of two zero- modes in super 
Hamiltonian of = 4 SUSY QM is caused by the fact that Witten theorem [22] modifies 
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when entanglement conditions are nonlinear (unlike the case of = 2 SUSY QM) as 
discussed in detail in [TS] . 

Let's consider the case when expression (|5]) contains one of the particular solutions, 
e.g. i^i{x,e), instead of ip{x,e,c). If one of the particular solutions of second order 
equation is known, the second solution could be calculated from the relation ip2{x,e) = 



ipi{x,e) J dt^- -. Thus, superpotential ([5]) obtains the form: 

J ipi[t,e) 



X 

W{x,e,X) = -hn{l + X J 



dt 



, (9) 



2 \(pi{x,e) + X(p2{x,e) I 2 \(p{x,e,X) 



It is easy to show, that the state with energy ii^ = in the spectrum of Hz is absent (i.e., 
the wave function tpz{x,E = 0) is nonnormalizable) and thus the spontaneously broken 
super symmetry exists. On the one hand, the state with zero energy and wave function 

ipl(x, £■ = 0) ~ — -Y, which could be normalized at certain values of A > 0, appears 

<^{x,e,X) 

in the spectrum of H^. In this case the exact N = 2 supersymmetry takes place. From 
the other hand, it is known [221 123] , that partial supersymmetry breaking is impossible in 
N = 4 SUSY QM without introduction of central charges. This contradiction resolves, 
when taking into consideration that employment of factorization energy e in construction 
of isospectral Hamiltonians is the simplest way of incorporation of central charges in 
= 4 SUSY QM. The similar situation occurs in consideration of form-invariant 
potentials [23]. More complete and consistent consideration of = 4 SUSY QM with 
central charges is given in [25J, but this is outside the scope of the current paper. 



3 General properties of multi-well potentials 

We will start the consideration of properties of isospectral Hamiltanians, derived in pre- 
vious chapter without concretization of the from of initial Hamiltonian, from obtaining 
the helpful mathematical expressions. Let yi and 1/2 be two linear independent solutions 
of homogeneous second order equation, then the following condition is hold: 



I 



W{yi,y2} 



{Am{t) + A2y2{t)y 



jdt 



Al + Al 



' A2yi{x) - Aiy2{x) 



^Aiyi{x) + A2y2{x) 



' A2yi{xi) - Aiy2{xi)'' 
Aiyi{xi) + A2y2{,Xi) 



(10) 



where W{yi,y2} = 1/1^2 ~ is Wronskian, which is independent on x for the second 
order equation, reduced to canonical form, as Shrodinger equation, and thus could be 
passed from integral. The special from of this relation for the parabolic cylinder functions 
was obtained in [10]. This relation is very useful for calculation of the integrals in the 
above given expressions. 
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Let's consider the case of c = 1 in the general solution of auxiliary equation 
This allows to avoid the unnecessary awkwardness, but nevertheless, reveal the funda- 
mental features of Hamiltonians and wave functions of = 4 SUSY QM, using only 
general properties of the solutions of auxiliary equation when e < Eq. First of all, it is 
natural to put Xi = — oo in expression ([5]) of superpotential, because (p{x,e,c = 1) = 

N^^ 1 
— tends to in this limit and function with asymptotic — — ?■ 

when X — ?■ — oo always exists, if using particular solution (pi{x,e). 



N- 



-2 



dt 



Ar-2 

[A(x, e, c) — A(— oo, e, c)] 



{l + C^)W{ip^,ip2}' ' ' ' ' ' ' ' " (11) 

Aix,e,c) = '-^P$^^^^ 
^i{x,e) + CLp2{x,e) 

Using above considered relations, it is easy to obtain the expression for superpotential up 
to constant term when c = 1: 

1 . , o r dt 



W{x,e,\) = --In l + AA^-' J 



'OO 



i(pi{x,e) + (p2{x,e)f 



1 + A(g,A)(^2(a^,g) A (12) 

2 \ ipi{x,e) + ip2{x,e) J 

A{oo,e,c = 1) - A - (A + 1) A(-oo,£:,c = 1) 
A(oo,e,c = 1) + A - (A + 1) A(-oo,£:,c = 1) 



Let's note that A(±oo,£, A), which enter the A{e,X), are determined by the asymptotic 
behavior of solutions of auxiliary equation. Due to this, while in the case of Hz potential 
is determined by symmetric combination (fi{x,e) + (f2{x,B), in the case of - by 
asymmetric (pi{x, e) + A{e, A) ip2{x, e): 

HZ = H-- ^ln((^i(x,e) + ^2{x,e)) (13) 

H+ = H--^\n i^x, e) + A(£, A) s)) (14) 

In some sense, potentials Uz{x,e) and U^{x,e,X) are form-invariant [15j, i.e. po- 
tentials and wave functions transform to each other by changing of parameters and 
their spectra are identical and this holds independently of the choice of initial Hamil- 
tonian. As it would be shown in the next chapter, if potential in Hz is multi-well 
symmetrical potential, then for H^ it is asymmetrical. Moreover, the form could be 
changed when varying e, as well as A. Relations (fTTj) and (fT2|) are useful also for 



7 



derivation of exact form of wave functions ipz{x,E) and ip+i^'-^)- Thus, expression 
for wave functions ip^{x,E) results from similar expression for ipz{x,E) by the sub- 
stitution ipi{x,e) + ip2{x,e) — )■ {(pi{x,e) +A{e,\) ip2{x,e)). In particular, the normal- 
ization constant for iIj^{x,E = 0) could be easily obtained from corresponding expres- 
sion for ijjz{x,E = 0). To demonstrate the possibility of partial symmetry breaking in 

= 4 SUSY QM we need to show that wave function iplix, E = 0) — — is nor- 

(^(X,£,A) 

malizable. Value of normalization constant could be derived from fITT]) when c = A > 0. 
Later we will return to the calculation of this constant for concrete form of initial Hamil- 
tonian. The potential in H^, in turn, has the same form as corresponding potentials in 
the case of exact supersymmetry with substitution A{e, A) — ?■ A. 



4 Isospectral Hamiltonians with almost equidistant 
spectrum 

To construct concrete expressions for potentials and wave functions we will choose the 
Hamiltonian with harmonic oscillator (HO) potential as initial. Let's consider the solution 

of auxiliary equation for e < Eq = —, {h = m = 1): 

' +2(£-^)U(x,£) = (15) 



Introducing dimensionless variables ^ = \^2ujx we obtain the equation for (p{^,€), where 
e 

u 

(^ + (^+2-TV V^(e,--) = 0,. = -2+^- (16) 
a) Exact = 4 SUSY 

This equation has two linear independent solutions: Djj{\/2^)^ Djj{—\/2^) - parabolic 
cylinder functions. According to terminology, introduced above, we denote V9i(^,e) = 

1^,(720, ^2{i.e) = D,{-V20 and Wronskian W{<fi,<f2} = EZ], where r(-z/) 

- gamma-function. Following the procedure from previous chapter, the general solution 
of auxiliary equation is chosen in the form: 

(^(e, e, 1) = D,{V20 + D,{-V20 (17) 

As one can see from f|T71) ip{x,e, 1) is an even function of ^. To obtain the exact form of 
superpotential, the integral, that enter the definition of the W{x,e, A) and normalization 
constant N~'^ need to be calculated from flTTl) . f[T^ . Due to the symmetry of (p{^,e, 1), 
expression of the integral simplifies and obtain the form: 



1 + \N-^ 



dt XN- 



+ A(+cx),£, r 



8 



Hence, superpotential could be expressed in the form: 

W{i^ e, A} = In f + + + + , (19) 



A(+oo,e,l) r(-z/) 

Using the asymptotic value of the parabolic cylinder function, we have A(+oo, e, 1) = — 1 
and expression for superpotential become quite compact: 

From (l20l) follows, that value of parameter A is limited by the condition A > — 1. Using 
the expressions from second chapter and superpotential (120|) for Hz and ipz{x,E), we 
could derive the exact form of Hamiltonians, which form the super Hamiltonian, and 
corresponding wave functions {ip^{C,,Ei) - wave functions of HO): 

HZ=H--^\n {D,iV20 + D,{-V20) , 

^-^^'^^^ ^ ^m^) ^i^) (21) 



Analysis of the potential fl2T]) reveals that several local minima exist only when < e < |. 
With values of e close to right boundary the third local minimum appears (FigJT]) and 
depth of outside minima increases. It should be noted that in the terms of dimensionless 
variable ^ the only way to vary the form of potential is through varying e and A. In the 
case of natural units, additionally, the form of potential(in particular, positions of local 
minima)could be changed through variation of u. Connection between H^, ip^{^,En) 
and H^, En) correspondingly could be treated in the same manner. Using fl20l) . the 

following expressions for H^ and ip^i^C^jEn) could be obtained from fl2Tl) by substitution 
ip{^,€,l) ^ ^{^,e,X + l): 

H+ = H---f-\n (D,{V2i) + (1 + \)D,{-V2i) 



i^ti^,E = 0) 



w{^^{^,E,),^{^,e,X + l)} 



2{Ei-e) (^(e,£,A + l) 



(22) 



{d,{V20 + (A + 1)D,{-V20) ^i^^ e,\ + l)' 



_2(A + l)v/^ 



"A+l 



Ti-u) 
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u 

0.30 
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Figure 1: Form of the potential U_ e, A = 0). 



From this relations the restriction — 1 < A follows. As could be seen from FigjJl potential 
U+{^i£i^) has well pronounced asymmetry, which increases when A — )■ —1. The value 
e = 0.47 corresponds the region, where U^{^,e,X) has three local minima. It shloud be 
noted, that depth of central minima also increases with A — > — 1. Given potentials and 
corresponding wave functions could useful for research of resonant tunneling phenomenon 

b) Partial supersymmetry breaking 

The utilization of the particular solution of auxiliary equation for derivation of su- 
perpotential (E]) realizes the situation of partial supersymmetry breaking. As it was 
mentioned above, the spectrum of Hz does not contain states with E = 0, because their 
wave function is nonnormalizable. Potential Uz{^,e) is presented on Fig|3l At the same 
time, zero state appears in the spectrum of with wave function: 

Normalization constant A^;^^ could be calculated using (|T9|) and has the form A^^^ = 
XW{ipi,ip2} = YJ~^^~j^ Thus, the limitation is A > 0. Hamiltonian has the form = 

^ <f _ 

— —— In {(p{^, e, A). Potential (e = 0.47) is presented in Figj2] at A > 0. 

When e = — - (-D_i(-\/2^) = e~y^| (1 — $(0); where $(^) - error function), poten- 
tials Uz{x,e) corresponds the potential, given in (the only difference is additional 
constant term). 



2 



2) \ 2 ) v^(l-<|.(0) 



2 



(24) 
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Figure 2: Potential e = 0.47, A). 




Figure 3: Potential Uz{^,e){0 <e < 0.5) 
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This potential is single- well, as well as corresponding ?7^(^,£: = —1/2, A), because exis- 
tence of several local minima is possible only when < e < -. Using the form-invariance 
property for U^: 

Uti^^ ^ = A) = ^ - ^ In ((1 + A) - (1 - A)$(0) (25) 

Spectrum of Hamiltonian with potential f l25p contains, unlike fl24p . the state with E = 
and wave function of the form: 

'^.«.^ = ")=((l + A)-(l-A)4K)) 



5 Conclusions 

The presented paper is devoted to construction and study of general properties of Hamilto- 
nians with multi-well potentials in the framework of = 4 SU SY QM. The emphasis is 
made on the research of the distinctive features of these Hamiltonians without concretiza- 
tion of the form of initial Hamiltonian. The given relation for certain type of integrals, 
containing the fundamental solutions of the Shrodinger type equations, allows to show 
the form-invariance of the isospectral Hamiltonians with multi-well potentials, obtained 
in the framework of = 4 SUSY QM. I.e, with the identical spectra of Hz and H^, 
corresponding potentials could be obtained from each other through the substitution of 
parameters. Moreover, this relation allows to analytically calculate normalization con- 
stants of zero modes, using only asymptotic values of fundamental solutions. Procedure 
of construction of isospectral Hamiltonians in the framework of = 4 SUSY QM also 
allows the partial symmetry breaking. On the model of harmonic oscillator we derive the 
exact form of isospectral Hamiltonians with multi-well potentials and corresponding wave 
functions. The existence of the parameters allows to vary the form of potential in the 
wide range. This is especially important in research of the phenomena such as tunneling, 
which are sensitive to the structure of the multi-well potentials. Moreover, the derived 
expressions could serve as a basis for construction of interesting exactly solvable models 
of stochastic processes, Kramers problem, stochastic resonance, to mention only few. 

Author is grateful to M.Plyushchay for valuable advices and appreciate M.Konchatnij 
and G.Ivashkevych for help with paper preparation. 
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